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Abstract 

Solvable vertex models in statistical mechanics give rise to soliton cel- 
lular automata at q = in a ferromagnetic regime. By means of the 
crystal base theory we study a class of such automata associated with 
non-exceptional quantum affine algebras Uq(g n )- Let Bi be the crystal 
of the Uq(g„)-modu\e corresponding to the Z-fold symmetric fusion of the 
vector representation. For any crystal of the form B — Bi ± ® • • • ® Bi N , 
we prove that the combinatorial R matrix Bm ® B — > B ® Bm is factor- 
ized into a product of Weyl group operators in a certain domain if M is 
sufficiently large. It implies the factorization of certain transfer matrix at 
q — 0, hence the time evolution in the associated cellular automata. The 
result generalizes the ball-moving algorithm in the box-ball systems. 

Keywords : Quantum integrable systems, Quantum groups, Crystal 
bases, Soliton cellular automata. 

1 Introduction 

The box-ball systems&EI are important examples of soliton cellular automata. 
They are discrete dynamical systems whose time evolution is expressed as a 
certain motion of balls along the one dimensional arrai£_of boxes. Their inte- 
grability has been understood by the ultradiscretizatkmEa of classical integrable 
systems (soliton equations). In the recent worksErEm it was revealed that the 
box-ball systems may also be viewed as quantum integrable systems at q = 0. 
Here by quantum integrable systems we jxiean the ones whose integrability is 
guaranteed by the Yang-Baxter equation,El and q is the deformation parameter 
in the relevant quantum group. In fact, the box-ball systems are identified with 
a q — > limit of some two-dimensional solvable vertex models, where the role 



of time evolution is played by the action of a row trans 
example is the original Takahashi-Satsuma automaton. 



■ matrix. The simplest 
whose classical origin 
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is the discrete Lotka-Volterra equation^ and the quantum origin is the fusion 
six-vertex model. Here is an example of the automaton time evolution. 

• ■ ■ 1112211211111111 • • ■ 

• ■ ■ 1111122121111111 • • ■ 

• ■ ■ 1111111212211111 • • ■ 

• - ■ 1111111121122111 • • • 

One regards 1 as an empty box and 2 as a box containing a ball. At each time 
step one moves every ball once starting from the leftmost one. The rule is that 
the ball goes to the nearest right empty box. One easily finds that the sequence 
of I balls propagate stably to the right with velocity I unless it interacts with 
other balls. By regarding such patterns as (ultradiscrete) solitons, the above 
figure illustrates how the larger soliton overtakes the smaller one with a phase 
shift. 

In terms of the fusion six vertex model at q = 0, the above hgure corresponds 
to the configuration: 
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The left and right boundaries are to be understood as 1 or 111 everywhere. This 
is a configuration of the fusion six vertex model in which the quantum space is 
spin 1/2 (1 or 2) and the auxiliary space is spin 3/2 (111, 112, 122 or 222). At 
q = only some selected vertex configurations have non-zero Boltzmann weights 
and the transfer matrix yields a deterministic evolution of the spins on one row 
to another. The vertex configurations in the above figure are the non-zero ones, 
and form an example of the combinatorial (q = 0) R matrix, which will be a 
main subject in this paper. Let Tm denote the row transfer matrix at q = 
corresponding to the spin M/2 auxiliary space. The above example corresponds 
to T3. Actually it can be shown that the ball- moving algorithm coincides with 
the action of Tm with sufficiently large M. In the above example Tm = T3 
holds for any M > 3. 

The coincidence of the ultradiscrete limit of soliton equations and the q — > 
limit of vertex models is an interesting phenomenon in various respects. From 
a statistical mechanical point of view, it roughly means that in those solvable 
vertex models, the profile of low-lying excitations over the ferromagnetic ground 
state at q = admits an exact description by (ultradiscrete) soliton equations. 
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From a mathematical point of view it leads to a systematic gencxalizationQ by 
means of the quantum affine algebras and the crystal base theory S3 See section 
4 of the referenceB for several examples of the scattering of ultradiscrete solitons. 

Now we turn to a general setting, where the six vertex model is replaced with 
a solvable vertex model associated with the quantum affine algebra U' q {Q n ). In 

this paper we treat the non-exceptional case Q n — A„ , A^_ x , A^ , B„ , , D\ 

and D^jl-y. The box-ball systems correspond to! g n = The Boltzmann 

weights are trigonometric functions satisfying the Yang-Baxter equation. The 
row transfer matrix is specified by the auxiliary space Vm and the quantum 
space ■ ■ ■ ®Vi j <g> Vi j+1 <g> • • • . Here Vm denotes the M-fold symmetric fusion 
of the vector representation. We suppose a ferromagnetic boundary condition, 
namely, the spins in the distance \j\ ^> 1 are all equal to some prescribed ele- 
ment in Vi r At q — the transfer matrix yields a deterministic evolution of the 
spins on one row to another. 

To analyze such a situation we invoke the crystal base theory Let Bi 
be the crystal of V. It is a finite set listed in Appendix [a| endowed with the 
action of Kashiwara operators ej,/, : B; — > Bi U {0} for < i < n. Let 
Si [a] € Bi be the special element as in (g). The states in the automaton are the 
elements • • • <g> bj ® bj+i <8> • • • £ • • • ® B\ j ® B\ ® ■ ■ ■ obeying the boundary 
condition bj — ^[ot], \j\ S> 1. Here afe is specified in Table [n] with (fi"l|), and 
k G Z is a label of the boundary condition at our disposal. Let Tm denote 
the q = transfer matrix, or M-th time evolution in the automaton; Tm ■ 
■••®5 fj ® B; J+1 ® • • • -> ■ • ■ ® (g> B ij+1 ® • •• . The map T M (- • • <8> & 3 - ® 
6j+i <8> • • • ) = • • • ® 6j <8 <8> • • • is induced by the isomorphism of crystals 
R : Bm ® B —> B ® B M with B = ■■■®Bi j ® B lj+1 ® • • ■ according to @. We 
call i? the combinatorial R matrix. It is obtained by successive applications of 
the elementary ones Bm ® B;, B; j <g) Bm- The boundary condition matches 

the known properties like <5m[<i/c] <8> ^ [ofc] >— ► ^ [afc] <8> #M[<Jfc]- When M gets 
large, Tm stabilizes to a certain map, which we denote by T. In the box- 
ball system terminology (§„ = Affl case), this corresponds to the boxes with 
inhomogeneous capacities {lj} and the carrier with infinite capacity. 

The main result of this note is Theorem 0, which states that the isomorphism 
R : Bm <8> B ^> B ® Bm with B — Bi x ® ■ ■ ■ ® Bi N is expressed as 

R={cr B ® a)PS lk+d ■ ■ ■ S lk+2 S lk+1 

in the domain BM[a/c] ® ^r'-j Bm <8 B ( |l2| ) if M is sufficiently large. Here Si 
is the Weyl group operatorcJ (|j) acting on Bm <8> B. P(u ® u) = v it is the 
transposition, it and <tb are the operators corresponding to the Dynkin diagram 
automorphism described around Proposition []]. The above result on R reveals 
the factorization of the time evolution in the automaton: (Corollary |l4|) 

T = (TBOi k+d ■ ■ ■ Si k + 2 Si k + 1 , 

where Si here is the one acting on • • • <g) Bi j ® Bi. +1 ® ■ ■ ■ . See the explana- 
tion after Corollary WM. Such a decomposition is by no means evident from the 
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defining relation (j32|). Note that T is a translation in the sense that the prod- 
uct <Jri k+d ■ ■ ■ ri k+1 {ji is a simple reflection) is so in the extended affine Weyl 
group, if a is interpreted as the Diagram automorphism acting on the weight 
lattice. According to the factorization, one can consider a finer time evolution 
T m (|3]) that includes the original one as %+td = T* (i > 0). For §„ = An \ the 
change from T m (p) to T m+ i(p) agrees with the original ball- moving algorithm 
in the box-ball systems, where one touches only the balls with a fixed color. In 
particular when = 1, our Definition |ll| provides a representation theoret- 
ical interpretation of the earlier observation!] For Q n ^ An , the automaton 
corresponding to VZj = 1 with afc = 1 has been introduced previously.Q The 
formula ( |33| ) in principle provides a simple algorithm to compute the refined 
time evolution for general lj and in an analogous way to the ball-moving 
procedure for the A^P case. The data d G Z>i and ii ^ / are specified in 
Tabl&f|. Curiously, they have stemmed from the studyE3c3 of Demazure crys- 
talsJHI^wilL-he-iuteresting to investigate the present result in the light of the 
worksBBO.BBE2 



2 Crystals 

Let fl„ = A^\n > l),^-!^ > 3),A%>{n > 2),B, ( l 1) (n > 3),C& x) (n > 
2),D n 1) (n > 4) and D^hOn > 2). For each Q n and I G Z> x , the U'{a n ) crystal 

B, has been" constmctedHEH except for C« with I odd. As for <#\ fl, here is 
B 1 '' in the paper.0 The finite set B; and the actions of ej, /j : B/ — > B; U {0} 
for i G 7 = {0,1, nr _rt} (crystal structure) have been defined. We employ the 
same notation asElll3 and quote the set Bi in Appendix [A]. In particular for 
Q n ^= An^ , we will always assume the convention 

x a = Xi, x a = Xi, a = i if a = i, 1 < i < n. 

Let us recall some other notations and the tensor product rule. 

<Pi(b) = max{j I fib + 0}, Ei (b) = max{j | e{b + 0}. 

For two crystals B and B', the tensor product B ® B' is defined as the set 
B B' = {&i ® 62 I 61 G B, &2 G B'} with the actions of ii and fi specified by 

ei{b l ^h) = {f^ l^\z Ei{ t\ (1) 
{ bi® eM if ^(01) < £i{b 2 )., 

{ h ® fib 2 if <p»(&i) < Ei(o 2 ). 



Here <8> 6 and 6 ® are understood to be 0. Consequently one has 
fi{bi ® 62) = ¥>i(&2) + ~ £i( b 2))+, 

£»(&! ® &2) = £i(&l) + (^i(&a) - 



(3) 
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where the symbol (x) + is.defined by (x)+ = max(x, 0). The Weyl group operator 
S*, (i € 2) is defined byN 



9h _ J /T l(6) - e<(6) 6 if >e,(6), m 
if^(6)<^(6). (4) 

Sj satisfies the Coxeter relations.0 For two crystals B and B', we let P : 
B®B'^>B'®B denote the transposition P(u ®v) = v ®u. It is easy to check 

SjP = PS* for any i E I. (5) 

For two crystals of the form B = Bi 1 ® • • • <& P; N and P' = Sj/ (g> • • • ® -B;' / , 
the tensor products B ® B 1 and B 1 ® B are isomorphic, i.e., they have the 
same crystal structure. The_isomorphism R : B ® B 1 — > B' ® B is called 
the combinatorial R matrix.U (In this paper we do not consider the energy 
associated with R.) It is obtained by a successive application of the elementary 
ones Bi i ® Bi>. By. <g> B^. We will use the same symbols ej, fi,ei,tpi,Si,P 
and R irrespective of the crystals that they act. 

Let A, denote a fundamental weight and let-P c i = ZAj be the classical 

weight lattice. (See Section 3.1 of the paperO for a precise treatment.) We 
define a linear map a : P c \ — > P c i as in the rightmost column of Table ffl. It 
is a Dynkin diagram automorphism. When g n ^ Cn (resp. jL = ck ), it 
agrees with the one introduced after Corollary 4.6.3 of the paperllil with B = Bi 
(resp. B = B 2 i) for any I. We also let a act on the index set / by the rule 



Table I: The data a' and a 



Bn 


a' on P/ 


(7 on P c ; 


A (1) 


a — > a — 1 


A a -> A a _i 




1 


A <-> Ai 


A (2) 


id 


id 




1 


A <-> Ai 




id 


id 




1 <-> 1, n <-> n 


A «-> Ai,A„ <-> A„_i 


D {2) 


id 


id 



i' = (r(i) cr(A,) = A^. For An \ the letters a, a — 1 should be interpreted 
mod n + 1. 

We also introduce the data d S Z >0 , and the sequences id, . . . , i\ G I and 
ad, ■ ■ ■ , cto for each algebra as in Table y. Here a^'s are taken from the letters 
appearing in the description of the crystals P/ in Appendix |X| In the third and 
the fourth columns of Table H the symbols and {-j-} mean that cither the 

simultaneous upper choice or the simultaneous lower choice arc allowed. 
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Table II: The data d, ik and 



Qn 


d 


i d ,. . . ,ii 




A {i) 

■tX.fl 


n 




1,2,... ,n 






1,2,... ,n + l 


^2n-l 


2n-l 
2n 


n— 1, . . 
n—1, . 


1 3, 2, {^q} , 2, 3, 
.. ,2,1,0,1,2,.. 


. . ,n 

. ,n 


n, . 


■ • '3,2, {y},2, . . . ,n-l,n,n 
n, ... ,2,1,1,2,... ,n,n 


o(l) 
U n 

n m 


2n-l 
2n 


n—1, . . 
n—1, . 


1 3, 2, | x } , 2, 3, 

.. ,2,1,0,1,2,.. 


. . ,n 

. ,n 


n, . 


■ ■ ,3,2, {-j-},2, . . . ,n~l,n,n 
n, . . . , 2, 1, 1, 2, . . . , n, n 


1J n\-\ 


2n-2 
2n 


n, n — 2, . . 
n—1, . 


• ' {1,0};^; ■ ■ ■ 
.. ,2,1,0,1,2,.. 


n — 2, n 
. ,n 


n, n 


-1,... ,2, {i},2, . . . ,n-l,n 
n, ... ,2,1,1,2,... ,n,n 



Proposition 1. For any I 6 Z>o i/ie operator a 1 := Si x ■ ■ ■ Si d : Bi — > is a 

bijection having the properties: 

Moreover, the action of a' is explicitly given by the second column of Table [J 

The second column of Table [j] specifies the transformation of those let- 
ters labeling the elements of Bi. See Appendix For example in case, 
a'((xi, . . . ,x n+ i)) = (x2, ■ ■ ■ ,x n +i,xi) in terms of the coordinates. Similarly 

for A 2n-V ' ■ ■ i x n-,X n , ■ ■ ■ = (xi,X 2 , ■ ■ ■ ,X n ,X n , . . . ,X 2 ,X 1 ). The 

proposition implies a' Si = S a [i)u' , from which the alternative expression a' = 
Si k+1 ■ ■ ■ S^Sg—i/i^ ■ ■ ■ S a -i( ik ) is available for any < k < d due to Sf = id. 
We identify a with a', thereby extend the definition of its domain. Namely, we 
also let a act on Bi for any I via a — Si k+1 ■ ■ ■ S^S^-i^ • • ■ Sg—i/^y (We do 
not exhibit I.) The result is independent of < k < d ancLenjoys the properties 
afi = f a {i)0 and aSi = e a {i^a. Proposition [l] was knownB for some k. For the 
tensor product crystal B = B\ x ® ■ ■ ■ ® B\ N , we write a B = ct !8> • • • <S> er : B ^ B, 
where a on the right side acts on each component Bi j of the tensor product 
according to the above rule a = Si k+1 ■ ■ ■ Si d S a -i(i 1 ) ■ ■ ■ S a -i^ k y Obviously one 
has a B fi = f a (i)0- B and a B Ci = e a(€j a B on B, therefore 

<r B Si = S a (i)a B on B = B h ® • • • <g) Bi N . (6) 

The combinatorial R matrix R : B ® B' — > B' ® B satisfies 

R(o- B 8) (t b >) = {vb> ® o- B )R. (7) 

When acting on crystals, a without an index shall always act on a single crystal 
Bi with some I. 

For each a E {1, . . . ,n + l}(g n = A}t'), {1, . . . ,n,n, . .. ,1} (§ n ^ vll 1 '), we 

set 

ft[o] = (a^a = i, x a , = for «V a) G (8) 
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with the notation in Appendix |A|. Using the crystal structure explicitly one 
easily finds 

5i[a k ] = S ik (<y,[a fc _i]) = e™ ax (*i[a fc _i]) 1 < k < d, (9) 
C-"C« = *H for any uG 5/, (10) 
<Pi k (8i[ak-i]) = 0, 

for any g„, where e™ ax 6 = e^^6. In particular (||) implies £;[<z<i] = cr~ 1 (Si[ao]) 
due to Proposition [j]. Thus it is natural to extend the definition of ik G I and 
a k to all k S Z by 

ifc+d = cr _1 (i)fe), 5j[afc+d] = cr -1 {Si[a k }) . (11) 

This way of extending the index k of a k is independent of I, and (g) also persists 
for all k 6 Z. We have {a/c | fc G Z} = {1, ... ,n + 1} for g n = and 
{1, ... ,n,n, ... ,1} for §„ ^A^. 

In this paper we will concern some asymptotic domain of the crystal Bm 
when M gets large. For a G {a k k G Z} and M ^> 1, we introduce the 
"domain" B M [a] C 5m by 

a - A (1) ■ 

Sn — Sill 

B M [a] = {(ui, . . . ,u n+ x) G B M \u a ^>u b for any b G {1, . . . , n + 1} \ {a}}, 

i a -il a » \u b - u b \ for any 1 
b G {1, . . . , n}, b^=a,a ( ' 

(12) 



-B M [a] = S (u%, ... ,ui) e B M 



In the rest of the paper we will have assertions under the conditions like c <gi 
u<g> c' G B (8 B M [a] <g> B' with B and B' of the form S = B h <g) • • • <g> B^,, 
B' = B Vi ®---®B v / (N, N' > 0). The mathematically invalid "definition" Q 
should be understood that the associated assertions are valid on condition that 
M ~S> 1 and the inequalities in ( |l2| ) are satisfied. Thus for example it amounts 
to assuming the following: 

e ik > ip lh on B(g> B M [a k -i] <E> B' , ip ik > e ik on B M [a k ] ® B' , (13) 

Ei k (u (g) x) = e tk (u) for u®x G BmK] 8 B, (14) 

^i fc+1 (x ® u) = <Pi h+1 (u) for x ® u <E B ® Bjfef[ofc]) (15) 

If u® x ^ y ® v then it (g> a; G -BmH ®B^>y®v£B® Bm [a] ■ (16) 

In the above (O) can be checked by using the explicit formulaBS for ^ and 
£j. (|l4|) and (^5p follow directly from ( |l5| ) and (|^). By the weight reason ( [l6| ) 
is obvious. Moreover we may effectively treat as 

S* 4fc (B <g> B M [a k -i] ®B') = B® B M [a k ] ® (17) 

for any fc G Z. 
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3 Combinatorial R matrices 



Let B be any crystal of the form B = Bi 1 ® ■ •• ® Bi N . Our goal in this section 
is to prove 

Theorem 2. When M is sufficiently large, the combinatorial R matrix giving 
the isomorphism R : Bm ® B — > B ® Bjf zs expressed as 

R={a B ® a)PS lk+d ■ • ■ S ik+s S ih+1 

on Bj«[(it] ® -B C £?m ® -B /or k £ Z. 

Definition 3. Let £> be any crystal of the form B = Bi x ® ■ ■ ■ ® Bi N . Given 

u <8> x £ Bm[oo] ® -B and y ® v £ B ® Bm[«o] we set 

u |fc> 0a;l fc > = 5« ifc . . . S ix (u <g> x) £ B M [a k ] <8>B, 
y W ® v ( k \ = S tk ■ ■ ■ S h (y ® v) £ B®B M [a k ], 

for < k < d. 

It should be noted that v) 1 *' for example is not defined solely from u but 
only with the other element x. 

It suffices to show Theorem ||for k = 0. To see this, let u(E>x £ B[ao] <8)B, 1 < 
k < d and assume u <8> x i— ► (<tb (8) a)PSi d ■ ■ ■ Si t (u(E> x) under the isomorphism 
Bm <8> -B ~ -B Cg) -Bm with M sufficiently large. Multiplying ■ ■ ■ Si 1 on the 
both sides and using (^]), one gets 

® x^ A (a fl ® a)PS ff -i (lt) • • • S CT - 1(il) S itJ • • • S lk+1 (u^ ® zl fe >). (18) 



In view of (|llj) and (Jl7|) this proves 1 < k < d case. To see it for the other k, 
multiply (cr m (g> ag 1 ) on the left side and (cr^ a m ) on the right side of @ for 
an integer m and use @ and (0). 

Henceforth we shall concentrate on the k = case of Theorem || in the rest 
of this section. Suppose Bm[o-o] ( S)BBu^)x^y®)v£B®) Bm[o.o] under the 
isomorphism Bm ® -B ~ .B <g) -Bm with M sufficiently large. Then the assertion 
of Theorem |^ with fc = is equivalent to 

v = cjv} d \ (19) 
y = a B x^. (20) 

We shall prove these relations separately. In our approach, ( |l9|) can be verified 
directly for any choice B — Bi 1 ® ■ ■ ■ <S> Bi N . On the other hand, as for Jib]), we 
first deal with N — 1 case and derive N general case based on it. 
Let us first treat (|l9|). 

Definition 4. For any crystal B we introduce 

t '. B — > ^>o 

6 .->(*! (ft),... ,t d (6)) 
t fc (6) = ^»(eS!* 1 -gr x 6)- 
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Here i^s are those in Table 0. It is easy to calculate t explicitly for B = B\. 
For the element b = (xi, . . . ,x„+i) for (b = (xi, . . . ,xi) for Q n ^ A^), 
the result is summarized in 

Lemma 5. The map t : Bi — ► Z> has the form: 
t(b) = (x n ,... ,Xx)for Ag\ 



(x n ,... ,x 3 ,x 2 ,{_ ' >,x 2 ,... ,i n _i) for A K 2 ^_ V 
(x n , . . . ,Xl,X ,Xt, . . . ,X n -i )for A™, 

(2x n + x ,x n -i,... ,X2,\_ 1,Xl >,x 2 , ... ,x n -x) for B^\ 

L^i, Xi) 

{x n , . . . , X\ , Xq , X\ , . . . , X n — 1 

) for CW, 

{X n + X n - U X n _2, ■■■ ,%2, S _ f ,^2, ■ ■ ■ ,X n - 2 ,X n -i + X n ) for 

\Xi,Xi I 



(2x )/or I? 



(2) 
n+1- 



See Appendix ^ for the notation Xq in -^2n ' ^ cllicl X0 in -Ci+i - The upper 
and lower choices correspond to those in Table [n]. From Lemma we derive a 
useful fact. 

Proposition 6. The map t : Bi — > Z> is injective for any I £ Z>i. 

Lemma 7. Under Definition |3| one /ias 

£. i( > |fc) ) = ^( u l fc - 1 > ^z^ 1 )), l<fc<d. (21) 

Proof. Definition | tells ® x^ = S lk {v} k -^ ® X \ k -^). Since u\ k ~^ ® 
^l*- 1 ) g _B Jlf [a / t„i] ® _B, this acts as 

®x^ = £l k u^ k -V ® elx^, (22) 
q = e lk (u^) - <p ih (x^ ) - (ul*" 1 ) ) - £ijt (^l^ 1 ) )) + , 

Thus we have e 4fc (u |fc) ) = £i fe (w |fc_1> ) -gS Vi k (u ]k ~ 1} ®a; |fc_1> ). □ 
For integers p, q depending on M in general, we write p < g to mean p < q 
or p — q = M-independent for A/ 3> 1. 

Lemma 8. Let u £ Bm[o>o] C -Bm /or sufficiently large M . For (ci, . . . , c,j) € 
Z| 0J de/me w< fc > (0 < fc < d) by = e** fc(u< * l)) ~ Ch u (l'-i) (1 < k < d) and 
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u' ' = u. Suppose cj < tj(u) for all 1 < j < d. Then the following hold for 



3 r^j L J 

1 < k < d 



( p ik (u^)=t k (u), (23) 
e lk (^ k) ) = c k , (24) 
t k (au^)=c k . (25) 



Although depends on c/s, the right side of ( |23] ) is independent of them. 
Actually (|24|) is trivial. The other relations in the lemma can be verified with 
a direct calculation by using the crystal structure of Bm- Under Definition ||, 
Lemma H immediately leads to (1 < k < d) 

cp ik (u\ k -V)=t k (u), (26) 

<Pi k (v {k - 1] )=tk(v), (27) 
e lk (u^) =t k (auW). (28) 

The right sides of (^6|) and ( |27| ) are independent of x and y in Definition [| 
respectively. 



Proof of ^1<\). Suppose Bm[o,o] ®BBu(3xt-^y(g>v£B(3 Bm[o-o} under 
the isomorphism Bm ® B ~ B ® Bm with M sufficiently large. We employ the 
notation in Definition^. Applying Si kl ■ ■ ■ Si t to the both sides of u&x ^> y(£>v, 
one gets u\ k ~ 1%1 <g) x\ k ~^ ni j/ fc-:l l therefore (p ik {v} k ~^ ® x^ 1 ^) = 

¥i k (y {k ~ 11 ® « <fc_1 ')- But 

are valid for 1 < k < d, telling that t(v) = t(au w ). Thus @ follows from 
Proposition ||. □ 
Now we proceed to the proof of (EOh with the simple choice B = B\. 



Lemma 9. Suppose £?m[oo] <8> B\ 9 <5m[oo] SmI^o] under 

the isomorphism Bm <8> -B; ~ ® -Bm M sufficiently large. Then we have 
t{u) = t(z). 

Proof. Define <E> z^ by Definition || starting from u®z = u'°^®z'°^ = 
<^M[ao] ® z. From ( |l9| ) we already know that u = au^. Thus we have 

t k (u) = t k (auW) i } e lk {u^) S ^(J*- 1 ) ® z^) 

( B^(^ |fe - 1> ) + (^(« |fe - 1> )-^(^- 1> )) + 
( B^(^ |fe - 1> ) + (4(u)- £i a^- 1> )) + - 
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Note from the explicit forms in Lemma [| that tk(u) =-ty(5M[ct>a]) = 0, from 

which tk(u) — ifi k (z' fe_1 ^) follows. In view of V4fc (t*l fc - 1 >) ^ t k (u) = 0and ([if), 
it also follows that z'^ = e™ ax • • -e™ ax z. Therefore we conclude (z\ k ~^) = 
ife(z) from Definition |^. □ 

Lemma 10. Given y (E) v e Bi <3 Bm[o.o], s ^t 

=g- ax ...g- ax ( 2/ ® w ) 

/or < < d. TTien we have t(av^) = t(y). 
Proof. Since = e™ ax • • ■ e™ ax y, we have 

v ( k ) = ~^ ik - iy )-^ k (v' k - iy ) v (k-i) = ^M"- x) )-^(y) v (k-i)_ 

Applying (fH we obtain t(av^) = t(y). □ 



Proof of (2C ) for B = B\. Suppose Bm[cio\®Bi 9 w(g)a; i— » y®w G I?i<8>i?M[ao] 
under the isomorphism £?m eg) Bi ~ (g> 5m with M sufficiently large. We 
employ the notation in Definition ||. Applying (cr Cg cr)e™ ax • • • e™ ax to the both 
sides and using (fio|), one finds 

B M [a Q } <8 Bi 3 5 M [a ] ® cra; |d) ^ (cr ® cr)g™ ax • ■ • £™ ax (y <g> u) G 5/ <8 B M M- 

Setting 

(CT ® Cr)^ laX ■ • ■ eX aX (y ® «) = (cry (d) ) ® (ot^), 

we have z = crj/ d ) and u = av^ d \ Thus t(u) = t(crt>w). But we know 
= i(o"xl^) from Lemma |^ and t(av^) — t(y) from Lemma [Tc|. There- 
fore i(crx' rf ^) = i(y), compelling y — ax^ by Proposition ||. □ 
To show (p0|) for the general choice -B = B^ <g> ■ ■ ■ ® , we prepare 



Definition 11. For s, s' G Z>o, 6, 6' G Bi and i £ J, we let the vertex diagram 

& 



6' 

denote the relations 

&' = ef db) ~ s)+ b, s> = Vi {b) + (s- £,(&))+. 

Here I can be any positive integer but we do not exhibit it in the diagram. 
The color i G I is attached to the horizontal line. The diagram should not 
be confused with the one representing the combinatorial R matrixia Given i, 
(6', s') is uniquely fixed from (o, s). Thus for example the diagram 
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b 2 



b N 



si ■ 



S-2 



SN-1 



S N 



is uniquely determined from sq, i and 61 
implies 



6jy- From Definition [ll] it 



b[ ® • ■ • ® b' N = g p— — + (6l g, . . . 8 6jv ). (29 ) 

Having established Theorem || for £? = £>;, we already know that 

B M [ao] <8> Bi 3 u (g) b £i cr& |rf) <gs cru |d> eB,® %[a ] 

under the isomorphism Bm ® -B; ~ -B; (g> -Bm for M sufficiently large, where 
u l fe ) = S ik ■■■S h (u®b) (0< k< d). 

Proposition 12. Under the above stated setting, the following diagram holds. 

b 



h(u) 
t 2 (u) 

td{u) 



il 






fell) 


il 






6 |2) 


id 


\d- 



t 2 {<T U W) 
1) 



b \d) 



Proof. By Definition [ll] we are to check 

6 i fc) = g K(^- 1> )-^(«)) +6 i fe _ 1>) 

t k (ouM) = ^ k {b ]k - 1] ) + (**(«) - ^(ft 1 *" 1 *)) 



for 1 < fc < d. To show the former, set x — b and apply ( p6| ) in 5' appearing in 
221) . The left side of the latter reads 



#^(^)) + (^(^-i))- £ifc(& l^>)) + , 

which is equal to the right side owing to (^6|). □ 

Proo/ of @) for B = B h (g> ■ ■ ■ (g> B 1n . Given any x = 61 ® • • • ® b N e B and 

it e B M [ao], set s ,fc = t k (u) (1 < fc < d). Let ftp e B ; ,s j>t € Z> (1 < fc < 
d, 1 < j < N) be the ones uniquely determined from the diagram: 
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b N 



SQ,1 
SO, 2 ' 

so,, 



ii 




in 








id 


b[ d - 



Si 2 



Sid 



ii 




ii 








id 





S2,l — 



S2,2 — 



S2,d — 



— SN-l,d 



ii 




1,1 












1,2 










b N 



SNA 



■ SN,d 



>N 



By a repeated use of Proposition fL2[ one has 
u ® 6i • • • ® b N 



ab^ <& cr&^ 



u 



under the isomorphism Bm ® B 
the other hand we introduce u' fc ^ 



-B <8> Bm with sufficiently large M. 



' fe > := S, 



■ 5i x (u x) G i?M[afe] 



On 

5 B. 



(Although not necessary here, v £ Bm[o>o] is uniquely determined from t{v) = 
(sjv,i ; ■ ■ ■ i s N,d) by Proposition ||, and we already know that the result coincides 

with v = en 
x \h) 



\d) 



from jis|).) We are to show 



,|d> 



b N can be proved by induction on k as follows. (We set 
<B> 6jv-) It is obvious for = 0. From (29) we have 



b' N ' . In fact, 



'N 



■Jk) 



'N 



_(m-s ,fc)+/i |fc-l) 



JV 



(30) 



where m 



\k-i) 



,\k-i) 

J N 



). On the other hand, x' k ' is determined 



from the recursion relation: 



+ Jfc-1) 



(31) 



because of £j fc (ul fe_1 >) » 1. Note that ip ik (u^ k ^^) 
two recursion relations ([30|) and ( |3l| ) lead to x^ = i 
induction assumption = 5[ fe_1 ^ ® • • • ® ^v"^- 

We have finished the proof of (|20| ) for any B = B[ 
the proof of Theorem ||. 



■Jk) 



tk(u) — so.fe- Thus the 



b N ^ under the 
□ 

N , and thereby 



Example 13. Consider q Ti 



3). The data in Table || reads ij 



cij = 4 — j mod 4 with ij G {0, 1, 2, 3} and a, G {1, 2, 3, 4}. To save the space 
the element [x\, X2,x%, X4) — (3, 2, 1, 0) G £?6 is denoted by 111223 for example. 
Then one has er(111223) = 112444 according to Table |. 

Let us take b <£) c = 111223 ® 344 G Bq <E> £?2_and seek its image under the 
isomorphism R : B 6 ®B 3 ^ B 3 ®B 6 . It is knownBthat R(b®c) = 223® 111344. 
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This can be derived by taking N = 1, M = 6 and k = 3 in Theorem p|, which 
reads R = (a (8 o^PS^-S^So. Namely we may regard 6®c€ -Bm[1] ® -B3, where 
M = 6 and xi = 3 for b are already sufficiently large so that 

b <g> c = 111223 ® 344 A 112234 ® 344 
A 112234^334 
A 112224 334 

^ 334 <g> 112224 ^5 223 <g) 111344 = i?(6 c). 

For comparison, take a smaller example i?(11223 <8 344) = 223 ® 11344 
B5. Theorem |^ is not applicable in this case under any choice of k because 
(11223® 344) = 11223<8344 for any i G {0,1, 2, 3} and (ct®ct)P(11223®344) = 
233<g)11244^ 223 ^ 11344. 



4 Cellular Automata 

The factorization of the combinatorial R matrix shown in Theorem |^ induces 
the factorization of the time evolution of the associated cellular automaton. 
Consider the isomorphism 

B M <g> (• • • ® B h ® B h+1 (g) • • • ) (• • • ® B tj <8 (g) • • • ) <g> B M 

induced by the successive application of combinatorial i? matrix Bm ® Bi j ~ 
Bj . ® Bm . We impose the boundary condition on bj G B^ as bj = Su [ak] for 
\j\ 3> 1, where the choice of k G Z is arbitrary. Assume the following properties: 

(i) Sm [ak] <8 ^[afc] 5i[a k ] ® ^ [a fc ] for any M, I, 

(ii) m (8> [a fc ] ® [afc] <8 • • • <8 ^. +JV [a fc ] ^ <8> • • • ® o i+JV ® flu* [a*,] 
for any u G Bm if N is sufficiently large, 

where bu (g) • • • <8>6; +JV is some element in i3; . <8 • • • <8 Bi +JV . ((i) is indeed valid by 
the weight reason.) Then under the isomorphism Bm®{- • ■®Bi j <8>Bi j+1 <E)- • • ) — 
(■■■®Bi. ®B l]+1 ®---)®B M , 

omM <8 (• • • (8> bj <g> ® • • • ) m- (• • • (gi bj <8 (8 • ■ • ) (8 <S M M (32) 

is valid for some 6^'s. One may regard the system as an automaton which 
undergoes the time evolution p = ■ • ■ (8 bj <g> 6 J+ i ® • • • 1— > p' = • • -®b'j (8 6^ +1 ® • • • . 
When M gets large, the transformation p 1— > p' stabilizes to a certain map, which 
we denoted by T(p) = p'. By taking £? = ••• <8 i?^. ® Bi j+1 ® • • • in Theorem || 
and using @, ( |Tl| ) we obtain 

Corollary 14. Under the boundary condition bj = 5i j [ak] G -0;- /or |j| 3> 1, 
£/ie time evolution T acts on B = ■ ■ ■ (8 B^ (8 -Bz i+1 (8 ■ • ■ as 

jii _ / a B^ik+td ' ' ' ^ik+a^-ffc+i if t ^ ^>0; 
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Actually all the Weyl group operators Si m for t > (resp. t < 0) in the 
above act as Si m = ejp ax (resp. Si m — /™ ax ) on B since they always hit such 
states p G B that e, m (p) 3> 1 (resp. </3j m (p) » 1). Corollary [1^ exhibits a factor- 
ization of the time evolution of the automaton having the background (vacuum) 
configuration •••(g) 5i j [afc] ® Si j+1 [afc] ® • • • specified by k e Z. When the partial 
factor Si m Si tn _ 1 ■ ■ ■ Si k+1 (k + d > m > fe + 1) in T is applied, the background, 
hence the boundary condition, changes into ■ • • ® 5/ . [a m ] (g) <5; j+1 [a m ] g) • • ■ ac- 
cording to @. A generalization of T that does not change the background in 
every intermediate step is constructed as follows: 



@k,m 



••Si 
•S 



fc + 2 + 1 ? 



m > fc + 1 , 



(33) 



where <7fc iTO acts on each component B\ i of the tensor product. We understand 
Tk = id. For any m > k, the operator T m retains the background in the orig- 
inal form ■ • • <g> Si d [an,] (g #i i+1 [fflk] (g • • • due to @. Moreover from Corollary 
|l4| we find that the original evolution under i-time step T* is included in {T m } 
as T* = Tk+td (t > 0). For A^P the evolution of the state p according to 
p = 7fc(p),7fe + i(p),7fc + 2(p), . . . agLfiea-ssdth that obtained by the ball- moving 
algorithm in the box-ball systemaHirEilEZI under a convention adjustment. In 
particular for , the evolution rule in terms of the Dyck language@ essen- 
tially agrees with the crystal theoretic signature rule in_applying Si — ef ax or 
ymax ^ _ explained in section 1.3 of the referenceJia 

Example 15. Consider §„ = A 5 (d = 5). To save the space the element 
ixi,x 2 ,x 3 ,X3,x 2 ,Xi) = (2, 1, 2, 1, 0, 1) G Bj is denoted by 1123331 for example. 
Let us take k = 0, hence the initial background is • • -igx^ [ 3 ] ®5i j+1 [ 3 ]<&••• . We 
employ is, . . . ,i\ = 2, 0, 1, 2, 3 and as, . . . , ao = 3, 2, 1, 2, 3, 3 correspondingly. 
Take 



P 



•33-12-3-31-2-33-33- - G ■ ■ ■®B 2 ®B 2 ®B l ®B 2 ®B 1 ®B 2 ®B 2 



where • stands for g) and ■ ■ • parts on the both ends represent the configuration 
identical with the background. Then the time evolution (|3l) is given by 



T{p) = ■ ■ ■ 33 • 33 • 1 • 32 • 1 • 23 • 33 • • • . 
According to Corollary O, the evolution is decomposed into the following steps. 



p = ■ 


•33 


12 


3 


3T 


2 


33 


33 


S 3 (p) = ■ 


•33 


12 


3 


31 


2 


33 


33 


S 2 S 3 ip) = ■ 


• -22 


13 


2 


3l 


2 


33 


22 


SiS 2 S 3 ip) = ■ 


• • 11 


13 


2 


32 


1 


33 


11 


S SiS 2 S 3 ip) = ■ 


••22 


32 


T 


32 


1 


33 


22 


S 2 S S 1 S 2 S 3 ip) = ■ 


•33 


33 


T 


32 


1 


23 


33 


a B S 2 S SiS 2 S 3 ip) = ■ 


•33 


33 


1 


32 


T 


23 


33 



Tip), 
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where in the last step we used the fact that ob interchanges the letters 1 and 1 
in each component as specified in Table |. Here the background configuration 
is changing except the last step. 

Alternatively the evolution may also be decomposed according to ( p3| ) into 
the following steps, in which the original background • • • ® 8u [ 3 ] (g> Si j+1 [ 3 ] (£> • • • 
is kept unchanged. 



p = %(p) = ■ 


•33 


12 


3 


31 


2 


33 


33 


Tx(j>) = - 


■33 


12 


3 


3T 


2 


33 


33 


%{p) = • 


•33 


12 


3 


2T 


3 


22 


33 


%{ P ) = ■ 


•33 


32 


1 


2T 


3 


22 


33 


%(p) = ■ 


■33 


32 


1 


32 


T 


22 


33 


%{P) = • 


•33 


33 


1 


32 


T 


23 


33 



= T(p). 



Regarding 3 as empty space in a box, one can interpret the above patterns 
as a motion of particles and anti-particles which can form a neutral (weight 
0) bound state. We hope to report on the explicit algorithm for general § n 
elsewhere. 
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A Parameterization of B\ 

We list the parameterization of the crystal Bi. In case, it may be identified 
with the set of semistandard Young tableaux of length I one row shape on letters 
{1,... ,rt + l}. For the other §„, Bi may be viewed as a similar set with some 
additional constraints. The relevant letters are {1, . . . , n,n, . . . ,1} as well as 
and/or depending on g n . The crystal structure (actions of e^, fi) can be found 
in the article! for Cn^ and the onesBEl for the other §„. 



n+l 



4 1} : Bt = {(. 



X\ , . . . 5 X n -\-\ 



)eZ n+1 \x t >o,J2^ = l}- 



n 




• B\ { (xi , . . . , X n , X n 



, . . . g z 2n \xi, Xi > o, yX x i + Xi ) = 



1G 



A 2n '■ B l = i( X l: ...,X n ,X n ,...,Xl) S Z 2 ™ | Xi,Xi > 0, y^(xj + x^) < 



i=l 



We set x — I — Y^(xi + Xj). 



Si 1 ) : = i ( Xl 



xi) e z 2n x {0,1} 



xo — or 1, Xj, Xj > 0, 



C$P : Bi = j(xi,...,x„,x„,...,xi) E Z 2 " 

n 

We set xo = (/ — Y^O^ + Xi))/2. 



1 > Efe=i(asi+^i) eZ-2Z 



i=l 



£>n ) : s i = i ■ ■ ■ ,x n ,x„, . . . ,xi) e 



x n = or x„ — 0, aci, Xj > 0, 
£" =1 (x;+x;) = 1 



-Di+i : S ; = \ (x 1 ,...,x n ,x ,x„,...,xi) eZ 2 ™ x {0,1} 



xq = or 1, Xj, Xj > 0, 

^0 + Yli=l( x i +Xi) <l 



We set xg = / — xo — T^(xj + i 
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